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Abstract: The Wishart distribution and its generalizations are among the most 
prominent probability distributions in multivariate statistical analysis, arising nat¬ 
urally in applied research and as a basis for theoretical models. In this paper, we 
generalize the Wishart distribution utilizing a different approach that leads to the 
Wishart generator distribution with the Wishart distribution as a special case. It is 
not restricted, however some special cases are exhibited. Important statistical char¬ 
acteristics of the Wishart generator distribution are derived from the matrix theory 
viewpoint. Estimation is also touched upon as a guide for further research from the 
classical approach as well as from the Bayesian paradigm. The paper is concluded by 
giving applications of two special cases of this distribution in calculating the product 
of beta functions and astronomy. 
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1 Introduction 

The Wishart distribution and its generalizations are among the most prominent probability 
distributions in multivariate statistical analysis, arising naturally in applied research and as a 
basis for theoretical models. The reader is referred to Gupta and Nagar (2000) and Anderson 
(2003) for a more extensive study regarding the theoretical as well as the practical uses of the 
Wishart distribution. Various generalizations and extensions are proposed for the Wishart dis¬ 
tribution, because of its importance in matrix theory. To mention a few: Sutradhar and Ali 
(1989) generalized the Wishart distribution for the vector variate elliptical models, however Teng 
et al. (1989) considered matrix variate elliptical models in their study. Wong and Wang (1995) 
defined the Laplace-Wishart distribution, while Letac and Massam (2001) defined the normal 
quasi-Wishart distribution. In the context of graphical models, Roverato (2002) defined the 
hyper-inverse Wishart and Wang and West (2009) extended the inverse Wishart distribution for 
using hyper-Markov properties (see Dawid and Lauritzen, 1993), while Bryc (2008) proposed the 
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compound Wishart and g-Wishart in graphical models. Abul-Magd (2009) proposed a gener¬ 
alization to Wishart-Laguerre ensembles. Adhikari (2010) generalized the Wishart distribution 
for probabilistic structural dynamics, and Diaz-Garcia and Gutierrez-Jaimez (2011) extended 
the Wishart distribution for real normed division algebras. Munilla and Gantet (2012) also 
formulated a special structure for the Wishart distribution to apply in modeling the maternal 
animal. 

There are of course many extensions that are not listed in the above, however the Wishart 
distribution can be viewed in the sense that it gives rise to other distributions. Thus the pos¬ 
sibility of extending each of the previous applications of the Wishart to hyper models, can be 
considered. We propose a possible construction methodology for creating new matrix variate 
distributions. The building block for our approach is discussed in the following section. To 
demonstrate the novelty, we compare it to a recent contribution by Carlo-Lopera et al. (2014) 
in the literature. 

Building Block 

Following Teng et al. (1989), Caro-Lopera et al. (2014) recently proposed a generalized 
Wishart distribution (GWD) under the elliptical models. They nicely derived the non-central 
moments of the likelihood ratio statistic for testing the equality of two covariance matrices under 
elliptical models for the corresponding matrices. Indeed, they considered the quadratic form of 
a matrix elliptical variate for building their distributions. We refer to p. 539 of Anderson 
(2003) and Diaz-Garcia and Gutierrez-Jaimez (2011) for more details and extensions. To be 
more specific, we recall a random matrix Y E ig gaid to have matrix elliptically contoured 

distribution with location matrix M G column covariance matrix S G Sm and density 

generator g : R'*' ^ R"*", denoted by ~ EC{M, S, 5 (), if its density function has the form 


f{Y) 

= \Y\-^g[tr{Y - M)'E-^{Y - Mf] 

(1) 


OR 


f{Y) 

= dn,m\Y\-hHY-M)i:-\Y-Mf], 

(2) 


where dn,m is the normalizing constant. 

Caro-Lopera et al. (2014) used Eq. (1) to develop generalized Wishart, however we deem to 
consider Eq. (2) in our construction. The difference in the form of the density generator g{.), 
plays deterministic role in extending matrix variate distributions. In Eq. (1), the normalizing 
constant is included in the form of g{.), however, it is not the case for Eq. (2) and the density 
generator in the latter equation can be any Borel measurable function. Thus, considering the 
quadratic form A = Y'^Y, the GWD based on Eq. (1) depends on the elliptical distribution, 
whereas the GWD based on Eq. (2) is free of any restriction and can take any form. The GWD 
based on Eq. (2) is neglected in the literature. This family of distributions, is a rich family with 
many applications. We propose some of the special members and applications in this paper. 

We organize the paper as follows: In section 2 a construction proposition behind the Wishart 
generator distribution is discussed using elementary tools in matrix theory and some of special 
cases are proposed. Section 3 contains some of the important statistical characteristic of this 
distribution, while a short note is given in section 4 regarding estimation purposes. Further 
developments beyond the Wishart generator distribution are given in section 5 and an application 
of a special case in section 5 is given in section 6. We conclude our result in section 7 and section 
8 is devoted to some necessary tools from matrix algebra. 
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2 Wishart Generator Distribution 


In this section a new family of distributions namely the Wishart Generator Distribution (WGD) 
is defined and some special cases along with the definition of inverse WGD are given. In a 
nutshell, the new generator type distribution concludes from a special case of Lemma [15] for 
K = 0 (see Appendix), which is provided in below. 


Definition 1 A random matrix X E Sm is said to have the WGD with parameter S E Sm, 
degrees of freedom n> m and Borel measurable function h{-), h{-) ^ 1 (called shape generator), 
denoted by X ~WGm((^,n, h), if it has the following density function 


/(X) = A:„,^|X|-t|X 
where using Lemma [75l for k = 0, 

rm(t)7o(f 


k~^ = 

'^n,m 


r( 


nm > 
2 > 


7o 


-h{tr'S-^X) 


y 2 h{y)dy 


provided that the above integral exists. 


Remark 1 The shape generator in DefinitionUl should sometimes admit the Taylor’s series 
expansion as a regularity condition, which will be referenced where ever needed. 


The reason of naming the distribution in Definition [Tj as Wishart generator, is the following 
result. 


Remark 2 Setting h{x) = exp{—^) in Definition{I\ yields the Wishart distribution (Press, 1982, 
5.1.1). Referring back to the the building block in the Introduction, it is clear that the form of 
h{.) here is free of taking any normalizing constant, however Garlo-Lopera et al. (2014) took an 
specific choice of h{.) to fulfill a valid density function for the Wishart distribution. 

Now we list some special cases, obtained from considering different selections of h in Definition 
[TJ Not to be conservative, various combinations of hyper geometric, trigonometric, exponential 
and Bessel functions can be considered to propose a new matrix distribution followed by WGD. 
The only restriction that should be fulfilled, is the existence of 7o(-)) Te., 7o(f) < oo. Looking in 
this way to construct a matrix distribution is not worthwhile from practical viewpoint, because 
it results in a complex structure. However, some applications are provided in section 6 for some 
special cases to address the practical importance. 


1. Taking h{x) = (1 + x) (2 +p) in Definition [T] for p > 0 we get the density function of a 
matrix variate t (MT) distribution as 


T ( r) T p) 1^1 I ^ m+1 , 1 . —( 

f{X) = ^ ^/J X|~2 |X|2—^(l + trX-^X) 2 


r,n(f)r(p)' 


( 3 ) 


where we used 


n 

70 [ 2 ) = 


nm 1 ^ t-. / nUl \ 

y 2 ^(1 + p) { 2 +P)dy = B \^—,pj 


2. Taking h{x) = exp (—ax^) in Definition [T] for a,b > 0, and using Eq. 3.478(1), p. 370 of 
Gradshteyn and Ryzhik (2007) we get the density function of a power Wishart distribution 
as 


f{X) 


ba 26 


T 


( nm \ 
\ 2b ) 


Isr2 |x|2 


m+1 

2 exp 


-a(trS-^X)'’ 


( 4 ) 
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3. Taking h{x) = (a + x) 2 exp(— 6 x) in Definition [H for \a\ < vr, 6 > 0, and using Eq. 
3.383(6), p. 348 of Gradshteyn and Ryzhik (2007) we get the density function of a matrix 
variate Kummer-type distribution as 


f{X) = 


\ — nm !— 

2^-Vb 




|X| 2 ' 


m+l 

2 


_ mn — L 

(a + trS~^X) 2 etr:{-b'E~^X) . (5) 


4. Taking h{x) = exp(— 6 x) (1 — exp(— 6 x))~^ in Definition [H for a < 1, 5 > 0, and using Eq. 
3.423(3^), p. 358 of Gradshteyn and Ryzhik (2007) we get the density function of a matrix 
variate logistic-type distribution as 

nm 

f{X) = ^ etr (- 6 X-IX) (1 - etr {-hY.-^X))-^ , ( 6 ) 

C7 ( 2 j 


where c = ^ • 

5. Taking h{x) = exp (—ax^) sin( 6 x) in Definition[Tl for a > 0, and using Eq. 3.952(7), p. 503 
of Gradshteyn and Ryzhik (2007) we get the density function of a sin-Wishart distribution 
as 


f{X) = 


nm+2 h_ 

2a 4 e 4 a 


bT iFi 


1 _ nm . 3. ^ 

4 ’ 2 ’ 4a 


X|t-^exp [-a(trX-iX)2] sm{bti'S-^X) , (7) 


6 . Taking /i(x) = exp (—x) ln(x) in Definition[Tl and using Eq. 4.352(4), p. 574 of Gradshteyn 
and Ryzhik (2007) we get the density function of a logarithmic-Wishart distribution as 

f(X) = * (-S-ljt) In (trE-'X) . ( 8 ) 


7. Taking /i(x) = pFq{ai, ..., Op, 6 i, ..., 6 ^; ex) exp (—x) in Definition [H for p < q, and using 
Eq. 7.522(5), p. 814 of Gradshteyn and Ryzhik (2007) we get the density function of a 
hypergeometric Wishart distribution as 


f{X) 


r(^) p+iF,(^,ai,...,ap, 6 i,..., 6 ,;c) 

|X|t ^ pFg (ai,... ,ap, 6 i,... , 6 g;ctrS"^X) etr (-X”^X) . (9) 


Theorem 1 Let X ~ WGm{'^,n,h). Then, Y = X ^ has an inverted WGD, denoted as 
Y ^ IWGm{'F!i,n,h), with the density 


f{Y) 


r (^mn) 
70(5^) rm(^n) 


det(S)-i’^ det(l")-^’"-5(™+i)/j(tr S-^l^-i). 


Proof: The result follows by the fact that under the transformation Y = X the Jacobian is 
given by J(X ^ r) = det(r)-(”^+^). ■. 

The following result gives some extensions to the existing result in the literature regarding 
matrix variate gamma distribution. 


4 













Definition 2 A random matrix Z E Sm is said to have the matrix variate gamma generator 
distribution (GGD) with parameters a > {m — l)/2, /3 > 0, X) € Sm, and shape generator h, 
denoted by Z ^ GGmi^, fd, h) if it has the following density 

7o(«)rm(a) 

Further ifW = Z~^, then W has inverted GGD with the density 

f{W) = det{'S)-^det{W)-^-^^^+^^h{2/3tv'S-^W-^). 

-fo[a)Tm{a) 

It is then denoted by W ^ IGGm{Sl,Oi,j3,h). 

Remark 3 Taking h{x) = exp {—^x) in Definition\^ gives the matrix variate gamma distri¬ 
bution of Lukacs and Laha (1964) and inverted matrix variate gamma of Iranmanesh et al. 
(2013). 

Note that if we take a = n/2 and (3 = 2, the GGD reduces to WGD. 


3 Properties 


Since the focus of this paper is the WGD, thus in this section we only give some important 
statistical properties of the WGD. These results can be directly derived for the IWD and GGD. 

It can be directly obtained that if W ~ WGmi'Sd, n, h), then the r-th moment of determinant 
of X is equal to 


E [det(X)n = 


mn) det(5])" 


Tm {\n) 7o {\n)^ 

r {j^mn) det(I])“5 


L 


det(X)"+5’^-5(”*+i)h(tr(5]-^X))dX 

n 

(r + \n) 7o (r + \n) det(S)'’+^"' 


r((r + in)m) 

r {\mn) Tm [r + ^n) 7o (r + ^n) 


T{[r + \n)m)Tm{\n) 7o 


■det(S)' 


( 10 ) 


Withers and Nadarajah (2010), demonstrated that for any square non-singular matrix X, the 
identity logdet(X) = trlog(X) occurs. Since logdet(X)'’ = rlogdet(X) = rtiTog(X) = 
trlogdet(X)'', from (fTO]) . for X ~ WGm(E,n,h) we have 


Eti log [det(X)^] 


^;iog [det(X)n 

log.B [det(X)n 

( T (Imn) Tm {r + ^n) 7o (^ + \ 

°^l^r((r + in)m)r^(in) 7o ) 


r log det(Xl). 


( 11 ) 


And if we put S = Im, then it yields 


Eti log [det(X)'’] 


.Blog [det(X)n 
log.B[det(Xf] 

( T (Imn) Tm {r + In) 7o (^ + !») \ 

°^\T{{r + ^n)m)Tm{^n) 7o / 
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Further for the expectation of zonal polynomial we have 


Emx)] = 


r (^nm) det(X)) 2 ’ 


'h) 7o ih^) 


det(X)5’"-5("*+i)C^(X)/i(tr(5]-iX))dX 


r (^nm) 


Tm (^n) 7o (^n) 


a (51). 


( 12 ) 


One of the important statistical characteristics of a distribution, might be its characteristic 
function (c.f). In the following result we give a closed expression for the c.f of the WGD. 

Theorem 2 Suppose that X ~iyGm(S, n,/i) and h admits Taylor’s series expansion based on 
zonal polynomials. The c.f is given by 

MT) ■ 

kT {^nm + k) 70 ( 2 ) 

Proof. The characteristic function is defined as 
V^^(T) = F;[etr(iTX)] 


I 


X VI 1 n , , n m+1 

^n,metr (zTX) |X| 2 |X| 2 2 


/i(trX-^X)dX 


— |X| 2 


etr (iTX) |X 


n _ m+1 


/i(trX-^X)dX 


(13) 


Using the Taylor’s series expansion we have 

tr(iTX)^ 


etr(iTX) = 

Hence, using Lemma [15] we get 


k=0 


k\ 


EE 

k=0 K 


Cn(iTX) 

k\ 


jL p 

MT) = EE ^pdet(S)-t / |X|t-'^C^(iTX)/i(trS-^X)dX 

7* n ... * Sm 


E E det(E)-tMdUil^ det(S)* C, (iTS) 

k=0 K 

00 -n / nm ^ 


r ( 7 ! + fc) 


EE 

k=0 K. 


r 


_lllLliM 

k%{i) r{!!f + k) 


C, (iTS). 


It might be ambiguous that how one can get the c.f of the Wishart distribution using the 
result of Theorem [2] Before rectifying this inconvenience, we need the following lemma which 
plays a key role in deducing the c.f of the Wishart distribution from Theorem |21 

Lemma 3 Let Y ~ IFm(El,n) (Wishart distribution of dimension m with n degrees of freedom). 
Then its c.f is given by 




k=0 K 
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Proof. It is easy to see that the c.f of Y has the expression iI)y{T) = det(/m ~ 2zTS)“'2. 
However in deriving the c.f, we make use of an integral over symmetric positive definite matrices. 
This integral is equal to 

Cl = j det(r)t-^etr +dl" 

= 2—det(5])2r™(^-jdet(I^-2iT5])-2, (14) 

On the other hand, writing the exponential term etr(zT’P) in CX as series of zonal polynomials, 
by Taylor’s series expansion, and using Lemma [15] for h{x) = exp (—we have that 


CX = [ det{Y)^-^ eti C^iiTY)dY 

k = 0 K V / 

OO - 

= EsE2”?‘*‘(f)/-(f)<iet(S)»C.(irE) (15) 

k=0 ' K ^ 

Comparing equations (fTTll and (IT^ . yields 

OO 

det(/m — 2iT5])~2 = 

fc=0 K 


which completes the proof. ■ 


Remark 4 Using TemmaQ it can be directly followed that by taking h{x) = exp 
remO we obtain the characteristic function of Wishart distribution, since ■'yj^ (^) 
and 7o(f) =2'^r(^). 


(—^x) in Theo- 

= 2=f+"r(^ + A)) 


Remark 5 If one is interested in deriving the distribution of a trace of a matrix, it can he done 
through inverting the Laplace transform. Using Theorem\^ replacing i by = —1, it can he 
directly deduced that the Laplace transform ofWG{'Y,n,h) is given by 


L{s) = Li[etr(—sX)] 

r(^)det(s)-t 


— li OO 


_/ nm 

S V 2 


+») 


70® 


k\ 


7o (i) det(sS )2 

r(^) 


EE 

A:=0 K 

^det [im - 


■c, (s-^) 


det (sE — Im) det(sS)' 


n+1 


7o{§) 

where the third equality obtained from Eq. (USD 


Theorem 4 Let X ^ WGm(Y,n,h), and A ^ Sm- ThenAXA'^WGm{{A') ^YA ^,n,h^. 

Proof: The proof follows from the fact that the Jacobian of the transformation Y = AX A' is 
given by J{X —)• Y) = det(A)“(™'+^). ■ 

It can be then concluded that if X ~ WGm{Y, n, h) and S = AA' , then AX A' ~ WGm{Irm n, h). 


Theorem 5 Let X ~ WGm {'Yi,n,h). The joint density function of the eigenvalues A = 
diag{Xi,Xm), Xm > ... > Ai > 0 o/X is given by 

“ /jW (0) 

k'.cAW 

det(A)t-'^A(A)C,(A), 
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where A(A) is the repulsion factor given by A(A) = A(Ai,..., Am) = 0 “ ^j)- 

i<j 


Proof. From Theorem 3.2.17. from Muirhead (2005), the density of A is given by 

— mP' p 

9 (A) = n / /(HAH')dH 

[t) r<j 


Note that from Dehnition[T] 


r( 


nm > 
2 > 


/ /(HAH')dH = 

Joim) 1 ml 2^01 2 i 

Since h{ ) admits the Taylor expansion, 


_n 71 m+1 

I] 2 A 2 2 


/ 0(m) 


/i(trS~^HAH')dH 


[ h(trS"iHAH')dH = f (S^IRAH') dH 

Jo(m) Jo[m) 


' 0(m) 

From dm, (Hi) and CSl) 

5(A) = 


fc=i K 


EE 

fc = l K 


h(^) (0) C, (E-1) a (A) 

(Im) 


—n (A. - A,) |S|-t 

Fm ( 2 j Tml 2 jTol 2 ) i<^j 

- /jW (0) (S-1) (A) 

X|A|2 2 > > — 


k=l K. 


k'.C^ (Im) 


( 17 ) 


(18) 


(19) 


„ m+1 ™ ™+l 

Note that [Al'z ^ = 0 ^ > hence 

2=1 


5(A) = 


7ri"*'r(^) 


r (^ 

^ m \ 2 




i<j 


)r™(§)7o(§) 

(0) C. (S-i) a (A) 


riY ’ EE 


2=1 


k=l 


klC, {In 


Theorem 6 Let X ~ WGm (X, n, h). Then for any A G Sm 

^{^) , , n hW (0)rm (?,k) Fm (^,^t) /ill 

P X < A =- .7,7 :sr 2 V V- V y m V 2 > ^^12 ^ / S-1A2 A 2 

^ ^ rm(f) 7 o(f)' ' fe!rm(§ + ^,«) "V 


Proof. Note that 

P(X < A) = 


r( 


nm > 
2 > 


. . TL t . . n 771+1 , 1 

r i-w / |X|2-—h(irS-iX)dX 

J^ml2j7ol2l 5 o<x<a 















Now make the transformation Y = A 2 XA 2 ^ then the Jacobian is J(X —Y) = lAI™'^^ 
hence 


P(X < A) = 


r( 


nm '' 
2 > 


r,n(i)7o(f) 

P(^) 


L 


|Y|2—^/i(trS~^A2YA2)dY 


0<Y<I„ 

( 0 ) 


r„^(§)7o(f; 


^ (0) f 


|Y|2- 


m+1 

■—a 


r( 


nm ') 
2 ) 


k=l n 
00 


0<Y<I„ 


S"^A2 YA2 dY 


rm(f)7o(t; 


y~! y~! c fs'^A^A^ 


fc=l re 


fc!r^(f+ ^,k) 


Remark 6 Note that X(^m) < a, is equivalent to 'K < aim since HAH^= X. To obtain the 
cumulative distribution function of \(^m)i largest eigenvalue of X. the previous theorem can 
therefore be used with A = aim- 


Theorem 7 LetX^ WGm{'X,n,h). The cumulative distribution function of X^m)^ the largest 
eigenvalue ofX is 




r(^) 


r,n(i)7o(t) 


EE 


/jl‘)(0)o“T„, 


(i-7 


■p ( m+1 

^m \ 2 ’ 




k=l 


k\Tm{^ + ^,n) 


Cn 


Proof. From Theorem 6 the cumulative distribution function of A(m)) the largest eigenvalue 
of X is 


^Nm) («) = 

T{^) 

rm(^)7o(7 

r(^ 


nri^i-EE 


k=l 

00 


hW (0) a^Tm{^,K) rm(^,Ar) 
fc!rm(§ + ^,K) 


Ck ( E ^0202 


_i-^i-a \ - \ - h^^^ (0) g^Fm (f, k) Fm k) 

n\I I 2-^2-^ 7.1V fn I m+1 n ) 


rm(f)70(t)' ' kWm{^ + ^,K) 

since (aS-^) = (X-^) . ■ 

Theorem 8 Let X ~ WGmiX, n, h). Then y = tr(X) has the following density function 
F(^)det(S)-texp(-y)^^ (f), 


fiv) = 


EE 


-It .,^+fe-i 


To(?) 

Proof. By applying inverse Laplace transform, using Remark [5l we get 


_ 1 \ riri 

S-i)y- 


1‘C+iT 

fiv) = / exp [sy] L(s)d, 

27rz r-)-oo 


fc =0 re 
00 fn) 

EE^e.(E-') 

fc =0 re 


7o(l) 

r(!fi) det(S)-5 

FTir 


rc+iT 

2vri 7c-iT 

1 


■ ( 2 +^) exp [sy] ds 


F (^ + A;) 


y 




( 20 ) 
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From Theorem[8l it can be easily concluded that for X ~ WGm{^,n,h), the r-th moment 
of tr(JC) has the form 




Ml) 


k=0 K, 


k^T{^ + k) 


( 21 ) 


In the following result, we derive the distribution of the ratios of the WGD in connection with 
the Wishart distribution. 

Theorem 9 Let X ~ WGm{o:'S,n, h) be independent ofY^ Wm{ldY,p). Then 
(i) The r.v. Bi = X-Wx- '2 has the following density function 

k r 

nn,m^ m \ 2 J I ^ \ 




2^r^(i) V/? 


det(S] 


p _ m+1 


EE 


k=0 

1 1 


k\ 




2(3 J r(m(^) + fe) 




(ii) The r.v. B 2 = {X + T)~ 2 X(X + T) 2 has the following density function 

k r (ti±P) /rv\ ^ 
m \ 2 J I ^ \ 


g2{B2) = 


2^r^(i) V/3. 

^ . n m + 1 ^ ^ p m+1 

det(B2) 2 2 det(/m-S2)2 2 

00 / V fc f n+p \ (n±P\ 

1 / a\ - Ir 


EE 

k=0 K. 


k\ 


2/3 J r(m(^) + fc) 


Proof: The joint density function of {X, Y) is given by 

f{X,Y) = Cdet(X)t-'^det(T)i-'^etr h/a-HiY-^X), 

where 


G = 




pm 

2^Tr 


(!) 


nm pm _ ^ n+p 

a 2 fd 2 det(S) 2 . 


Make the transformations Bi = X ^YX 2 and U = X with the Jacobian J{X,Y 
BuU) = det(I7)5(™+i) to get 


g{Bi,U) = f{U,U2BiU 


^ , n+p m+1 , , p m+1 / I il l\.i is 

= Gdet{U)^ -^ det(Si) 2 —^ etr fS-^t/ 2 Bi ?72 J h{a-^Y~^U). 

For a moment, assume that the distribution of Bi is symmetric. Thus from symmetrized density 
we have 


giiBi) = 


[ gi{HBiH')dH 

J 0(m) 

[ [ g(HBiH',U)dUdH 

J@(m) Jsm 


, s P m + 1 / , , n+p m+1 1 , 

= Cdet(Bi) 2 —^ / det([/)^- ^h{a-^Y-^U) 

Sm 
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[ etil-^U^'E-^U^HBiHAdHjdU 

J&lm) \ 2p J I 


I0{ 

= Cdet(B] 


p _ m + 1 


EE 


/ det{U) — 

By making use of Lemma fT5l we get 


= Cdet(Si)i-^ 


k=0 K 
n+P _ m+1 


k\ 


1 yc^iBi) 


2(3) C^{iy 


C^{T.-^U)h{a-^T,-^U)dU. 


EE 


fe=0 K 


fc! 




2/3 y C^{Im) 

^ ^ ^ { -det S ^C+I. 


r(m(^) + fc) 


(Im). 


After simplification, we obtain (i). For (ii), make the transformations B 2 = {X + 1^) 2 X(X + 


l")-2 and F = X + 1", with the Jacobian J{X, Y B 2 , V) = det(F) 
g2{B2) = Cdet(S2)2-^det(I^-B2)i-^ 


m+1 


to get 


L 

%J On 


. n + p m+1 

det(F) 2 2 


1 


etr ( -j^V^Y-W 2 [Im-B 2 ] ) /i(a“^ tr F 2 S"! V 2 B 2 )dF. 


m+1 


Make the transformation Z = Y 2 with the Jacobian J{V Z) = det(S) 2 and use 
Lemma [12] to obtain 

52(^2) = Cdet(B2)§-'^det(/™-S2)§-'^det(S)'^ 


^ , n+p m+1 

det(F)^- 


1 


= Cdet(S2) 

00 ^ 

EE^ 


p _ m + 1 


2 etr (--B 2 ] j/i(a ^tvZB 2 )dZ 

p m + 1 , n+p 


det(/m - B 2 


det(S)- 


k=0 K 


1 


{Z[Im - B 2 ]) h{a-^ tr ZB 2 )dZ 


f det(F) 
t(lm — B 2 ) 

“ _ (!!± 2 ), r„ (!!± 2 ) 7 . (!“) o”'('P)+‘ 

k!\ 213 


, ^ p m + 1 , , p m + 1 , , Ti+p 

= Cdet(S2)2—^det(/,^-S2)2—^det(S) — 


fc=0 K 


r(m(^) + A:) 


n + p 1 

det(B2)-—C,(S2 


After simplification, gives (ii) and the proof is complete. ■ 

Remark 7 One way of checking the accuracy of the result of Theorem\^ is to consider whether 
one can get the same result by taking h{x) = exp (— for the Wishart distribution. It is well 
established that if X Wm{'Zl,n), then Bi has the well-known beta type II distribution. This 
result directly follows by making use of Eq. (HSj). It can be also shown that B 2 has the beta type 
I distribution if we take h{.) to be of exponential form. 


4 Estimation 

In this section, we briefly consider some estimation aspects for the WG distribution, including 
the classical as well as Bayesian viewpoints. The focus is the latter paradigm. 
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4.1 Maximum likelihood estimation 


In this section, we derive a non-linear equation to find the maximum likelihood estimator (MLE) 
of X) along with Fisher information matrix. 

Theorem 10 Let X ~ WGm (X, n, h), where the trio {m, n, h) is assumed to be known. Further 
assume that h{-) is a monotonic continuous and differentiable function. Then the MLE ofS is 
given by 

± = ‘^g' (tr(5]”^X)) • X, 
where g{.) = — log[/i(.)] and g'{x) = dg{x)/dx. 

Proof. The likelihood function is given by 

L (S) oc det(X)-t det(X)t-"^/i(trS~^X) 

Hence the log-likelihood function is 

/ (S) a ^ log det(S-i) + log det(X) + log [hitr'S'^X)] 

cc ^ logdet(S~^) -|- log [/i(trS“^X)] 

To obtain the maximum of the log-likelihood function, let Z = trS”^X; then differentiated 
log-likelihood function has the form 

cx |[2X-diag(X)]-^^[2X-diagX]. 

Setting to zero gives the MLE of S as 

X = -g'{Z)-X 

n 

= ^9'(tr(E-'x)) - X, 

■ 

Since the structure discussed in Theorem IIUI is similar to the generalized elliptical distribu¬ 
tions studied by Frahm (2004), we do not provide inferential aspects of the MLE here and for 
complete explanations on the MLE regarding existence, consistency, applications and etc., the 
reader is referred to Frahm (2004). 


4.2 Bayesian estimation 

Theorem 11 LetX|S WGm {'S,n,h). Suppose that the prior distribution ofT, is an inverse 
Wishart distribution with parameters FI andp, hence S ~ Wff^ {Fl,p). The marginal distribution 
of X is given by 


m(X) = 


r(^) 


p(p — m — 1 ) 


r^(i)r™(t)7o(f) 


^ ^ p —m—1 

det(fi) 2 


X det(X) 


m + 1 ^> 

fc=l K 


r (m^ + k) 


a. 


2 
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Proof. The marginal distribution of X is given by 


m(X) = [ /(X|S)7r(S)dS 

JSm 


r( 


nm > 
2 > 




, , n m+1 ^ , p — m — 1 

det(X) 2 2 det(fJ) 2 


det(S)- 


n _P 
' 2 2 


etr ( ) /i(trX"^X)dX. 


Now, let S ^ = T then the Jacobian is det(T) 


—m—1 


m(X) = 


r( 


nm > 
2 > 


2"^^r™(|)r^(f)7o(t) 


hence 


, , n m+1 ^ , p — m — 1 

det(X) 2 2 det(fJ) 2 


f det(T) 

Sm 


r( 


n+p —m— 1 
2 


nm > 

2 ^ 


2"^^r™(|)r^(f)7o(t) 


etr ( /i(trTX)dT 

, , n m+1 ^ , p — m — 1 

det(X) 2 2 det(fJ) 2 


00 « 

det(T) 

k=l K 


n + p —m — 1 / I \ 

- 2 -( --Tfi j /i(trTX)dT. 


Using Lemma [15] we get 


det(T)^^^^V^C« (-^Tn) /i(trTX)dT 


Hence 


r (m^ + k) 
r{^) 


n+p 


det(X)-^C, ( 


p(p-m-l) 

2^^r^(i)r™(§)7o(§) 


^ ^ p — m—1 

det(ri) 2 


CO /'IL+Pt r tli+Pl tL+Pf 

xdet(X)-i-'^ E E ^ ^ ^ ( -7nx 


fc=l K 


r (m^ + fc) 


Theorem 12 Z/etX|S ~ WGm (X,n, h). Suppose that the prior distribution o/S is an inverse 
Wishart distribution with parameters SI andp, hence X ~ posterior distribution 

of S is (/ipen by 


7r(S|X) = det(X)®Pdet(E)“t“5etr ft(trE"'X) 


g ^ Liox-- 


U.1 « r(m7£ + i) 

Proof. The posterior distribution is from Bayes’ theorem as 

/(X|S)7r(S) 


n -1 


7r(X|X) = 


m(X) 
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Hence 


TrfSlX) = 


2"^^r^(|)r^(f)7o(f) 


, , n m+1 ___ p-m- 

det(X) 2 2 det(fi) 2 


xdet(S)-t-ietr f/i(tr5]-^X) 




00 ( n+p \ p ( n+p \ ( n+p \ 

I 2 2 )lk\ 2 ) 

r(m.7£ + t) 


y—HL — X , 'P \ m-h. 

det(n)-^ det(X) 2 +— 


c,(--nx- 


n+p n p / 1 1 \ 1 

= det(X)—det(X)-2-2 etr f h{tr^~^X) 


Jf-V r Z' 


Theorem 13 Lei X|S ~ WGm (X,n,/i). Suppose that the prior distribution ofT, is an inverse 
Wishart distribution with parameters S7 andp, hence X ~ {il,p). Then the Bayes estimator 

0 / IXI under the squared error loss function is 


1.1=1 l^x 2 “^ ) 


E oo 

k=l T. 


('n+jp\ r ! 'n+p \_, ( n±p\ 

V 2 )k^\ 2 /7fcV 2 / /-< f lOY —1 


det(X) 


r[mll±^+k) 




Proof. The Bayes estimator of det(X) under the squared error loss function is 


det(X) = L;[det(X)|X] 


J det(X)^det(X)-t-|+ietr /i(trX-^X) 

t^i. T{m^ + k) V 2 J 




00 « 

det(T) 

1=1 A 


n+p—2 —m —1 / I \ ^ 

- 2 - Cx —TQ /i(trTX)dT 


where T = X Note that 


f ^ ^ n+p — 2 — m — l 1 \ ^ 

/ det(T)- 2 - Cx ( “ 2 ^^ ) ^(trTX)dT 

( 72 -l)>r„(!!± 2 -l) 7 , / 1 

- L\x\-^+'cJ--nx-' 

r — m + 1) 72 ! 
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from Lemma [T5l Hence 


det(S) = det(X) 


OO 

EE 

.k=l K. 


( n+p \ -p 


/ 1 


r (m^ + k) 


a, I "2^^' 


-1 




1=1 A 


r (m^ -m + k) 


2^1=1 l^X rfmlliii-m+fc'l 2 “^ ) 


r(mll^-m+k) 

( n+p \ p / H+El^, 1 H+El 
I 2 )k'^\ 2 nk\ 2 ) i lr»v-l^ 
l^k=\2-.K r(m=^+fc) 2 “^ ) 


■det(X) 


5 Further Developments 


In this section we provide the reader with some plausible extensions of WG distribution. In this 
respect, we first define the hyper geometric WGD as in below. 

Definition 3 A random matrix X E Sm is said to have the hypergeometric WGD with param¬ 
eters ai,..., Op E C, bi,... ,bg E C, (p < q), D, S E Sm, degrees of freedom n > m and shape 
generator h{-),h{-) ^ 1, if it has the following density function 

/(X) = /„,mdet(X )"2 det(X)2 ^ pFq (oi,... ,ap;6i,.. .,bq;ftX) /i(tr S~^X) 


where 


/-I 

'"n,m 


det(X) 

det(X)-t EE 

fc=i 

[ det(X)t 


t / det(X)t-"^ pFg(ai,...,ap;6i,...,6g;DX)/i(trS-^X)dX 

^ Sm 

°° (ai)^,..., (ap)«, 1 


k=0 K 

n _ m+1 


JSm 

det(X)' 


= r. 


(^l)re) • • • ) (bq^K, k\ 

/i(trS-iX)C«(DX)dX 

(Qi)k) ■ ■ ■ ) (Qp)k 1 ('2)k ^"^(2 )7A:('2) 

{bl)^,...,{bq)^k\ T{r^+k) 

id E E I?!'" " ’ c-(ns) 


5EE 

OO 




nm 

2 


+ A:) 


det(S)2G«(DS) 


/rom Lemma 1731 ILe designate this by X ^ HWGm{'^,^, 0 ',h,n,h), where a = (ai,...,ap) 

and b = (bi,... ,bq). 

As a direct consequence of Definition [3l taking p = 0, q = 1, bi = ^, Q = for E Sm, 

gives the non-central WGD as in below. 

Definition 4 A random matrix X E Sm is said to have the non-central WGD with parameters 
G Srri) degrees of fveedoTu n ^ in and shcipe geueratov ^ denoted by X ^ 

A^VrGm(S, n,/z), if it has the following density function 

/(X) = ln,m det(S)-t det(X)t-'^ oi^i Qn; /i(tr X’^X) 
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where the normalizing eonstant is given by 


'"n,m 


= r. 


oo 

i)EE j 


k=0 K 


Ikin) 


kT{^ + k) 


CJ^) 


since 

Another interesting distribution raises from Definition [3l comes up by setting p and g to 0 and 
1, respectively as: 

/(X) =/„,^det(S)-?det(X)?-^etr(nX)/i(trX-^X). (22) 


We call this distribution as the exponentiated WG distribution. Note that according to Theorem 
[T^ the posterior distribution of S has the exponentiated WG distribution. 


6 Applications 


In this section, we briefly consider some applications of two special cases of WGD. 

Distributions of the form ([3|) has many applications. Arashi et al. (2013) showed that the 
posterior distribution of scale matrix in the matrix variate t-population under Jeffreys’ prior has 
the MT distribution given by (l3|). Another interesting application of the MT distribution is the 
following result, where we show that finite product of beta functions can be written as a ratio 
of gamma functions. 


Theorem 14 Let p > m, then 

n+l 


nB(-.p+(i-2)-)=h#hJh 

ii (2’-^ ^ ’2) T{^+p) 


i=0 


Proof. Using Gorollary 3.2.3 of Srivastava and Khatri (1979) for Z = (Zi,..., Z„), Z, G 
we have 


/ , ^ ri m+1 ^ ^ { nm , 

/ det(X) 2 —^ (l + trX)“l~+^’ldX 

I 

Jr^ V i=i 


-i^+p) 


Ui=Zj Zi 








diZ \... dZji 
-(2f+p) 


Ui 


dui . . . dUr, 


i=l 




Ui 


-i^+p) 


+ Ul 


dtti... dun- 


Now apply the transformation Vi = for i = 2,... ,n, with the Jacobian J(tt 2 , ■ ■ ■ ,Ur, 

V 2 , ■ ■ ■, Vn) = (1 + ni)"’“^ to obtain 


lb lb 

U2 = U2(l +Ui), ^ = (1 + \ 


i=2 


i=2 


Hence we get 


I = 


J{0,1) 


16 







X 


n 


--1 


^ n yi^+p) 

1 + ^ f i dV2 ■ ■ ■ dVr, 


i=2 

B(^,p-m 


i=2 


--1 


Un 


'( 0 , 1 )" 




-(=f+p) 


i=3 


i=3 


l + Vi 


dv2 ... dt>„ 


Again make the transformation Wi = for i = 3,...,n, with the Jacobian J{v 3 ,... ,Vn 

W 3 ,..., Wn) = (1 + V2Y~'^ to get 


1 



dW3 . . . dWn 


Continuing this procedure, finally yields 


1 


n+l 

ns 


i=0 



But since f{X)dX = 1, from Eq. ([SD for E! = I^n, we have 


I 


, ^ n m+1 , ^ f nm i 

det(X) 2 —^ (1+ trX)"i~+P-' dX 


rm(f)r(p) 

T{^+p) ’ 


which by substituting in Z, completes the proof. ■ 

For considering another application, let Y ~ EC{M,'Y,g) and consider the distribution 
oi Z = Y^Y. It is well-known that if Y has matrix variate normal distribution, then Z has 
Wishart distribution. For a moment let T = S “2 . Anderson and Fang (1982) 

derived the density of Z for the case M = 0 and T = I^- Fan (1984) extended their result by 
presenting the density of Z for general M and T as an integral form. Afterward, Teng et al. 
(1989) derived the closed form of the density Z for practical use. 

Now as an application, we show that the distribution of Z is the non-central WG. To see this, 
consider that using Theorem 1 of Teng et al. (1989), if "F ~ EC{M, 'Yi,g) then the distribution 
of Z = Y^Y is given by 


fiZ) = 


mn 

TT 2 

Tm (f) 


|S| 2 |2’| 2 


m + 1 
2 


k=0 


A2k) 


(tl 


(S-^Z + T)) 


k\ 


E 


& TE- 


1 1 

'2 ZYi 2 


(§), 


(23) 


where g^‘^^\.) is the 2A:-th derivative of g{.). 

If we take p = 0, g = 1, 6i = (|), $1 = S“ 2 TS “2 and h{x) = g^‘^^\x J-trT), then using 
Definitionll Z ~ lVlTGm(S, fi, n, h). 

For an application of the NWGm{Y,fl,n,h) was introduced here, consider the use of the 
non-central WG distribution when it arises from lighter/heavier marginal tail alternatives to the 
matrix variate Gaussian distribution, in astronomy (see Feigelson and Babu, 2012 for applica¬ 
tions of statistics in astronomy). To be more precise, in the study of imaging extrasolar planets 
for life, as discussed by Tourneret et al. (2005), direct imaging through statistical signal pro¬ 
cessing is the only method for exoplanet detection. See Figure 1 (adopted from Google Images) 
to set the platform for investigating the true distribution in the forthcoming explanation. 
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PLANETAS EXTRA-SOLARES 



Figure 1: Visualizing extra-solar planets from their position distribution, as it might be seen 
from telescope. 



Figure 2: Instantaneous intensity of the wave in the focal plane. 


As explained by Aime and Soummer (2004) the complex amplitude of a wave in the focal 
plane of a telescope, at a position (x,y), can be written as follows: 

= C(x,y) + S{x,y), 

where C{x,y) G C is a deterministic term proportional to the wave amplitude in absence of 
turbulence and S{x,y) G C is the wavefront amplitude (associated to the speckles) distributed 
according to a zero mean complex Gaussian distribution. Tourneret et al. (2005) assumed 
that the telescope aperture has central symmetries which imply C{x,y) G M and using the fact 
that the real and imaginary parts of 'ip{x,y), denoted by il)^{x^y) and ijji{x,y) have Gaussian 
distributions, extended the instantaneous intensity of the wave in the focal plane at a position 
ix,y), given by 

A(x, y) = |V’(x, ?/)p = V'^(x, yf -h 'il^iix, yf 
to multidimensional case. They demonstrated that 

A = (A(l), . . . , = V’rV’r + 

where -0^ = (V^^(l),..., V'^(n^))^ and V'r = (V’i(l)) ■ ■ ■ ? : has non-central Wishart distri¬ 

bution. 

Since speckles are bigger than they appear in telescope, it is highly misleading to assume the 
normality assumption, even if the assumption of symmetry is taken, to study of planet formation. 
See Figure 2 for the distribution of amplitude of wave in the focal plane of a telescope. 

Thus it is more plausible to take these speckles as extremes in amplitude of a wave or 
outlier in plane formation as appears in telescope. In conclusion, accepting the assumption of 
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symmetry, the multivariate t-distribution (or may be lighter tail alternative to Gaussian, as it 
might be captured from Figure 2) is a relevant alternative to the normal one. Hence, by the 
theory discussed in the above, A = (A(l),..., A(n^))^ has non-central GW distribution arises 
from taking g to be the kernel of multivariate t-distribution in Eq. (2). 

7 Conclusion 

In this paper a family of distributions were introduced from the Wishart generator distribution 
which includes the Wishart as a special case. The Wishart generator distribution might be 
important for a number of practical signal processing applications including synthetic aperture 
radar (SAR), multi-antenna wireless communications and direct imaging of extra-solar planets, 
the latter was discussed using the non-central Wishart generator distribution. Several statistical 
properties of this newly defined distribution were studied from matrix theory viewpoint. Brief 
notes regarding classical as well as Bayesian estimations were also proposed. 


8 Appendix 

Initially let Sm and Im be the spaces of all positive definite matrices of order m and all symmetric 
matrices between 0 and Im under the meaning of partial Lowner ordering, respectively. For a 
given matrix A G denotes the transpose of A, tr(A) = an app\ etr(A) = 

exp(tr(A)); det(A) = determinant of A; norm of A = ||A|| = maximum of absolute values of 
latent roots of the matrix A; and A 2 denotes the unique square root of A. 

Also denote the space of all orthogonal matrices of order m by 


0(m) = {H\H'H = HH' = /^} , [ dH = 1. 

J 0(m) 


Let A: be a positive integer; a partition k of /c is written as k = (ki, k 2 , ■■ ■), where kt = k. 


Definition 5 (Muirhead, 2005) Let Y be an mxm symmetric matrix with latent roots yi, • • • , ym 
and let K = {ki, ■ ■ ■ ,km) be a partition of k into not more than m parts. The zonal polynomial 
ofY corresponding to k, denoted by C^iY), is a symmetric, homogeneous polynomial of degree 
k in the latent roots yi, • • • ,ym such that: 

(i) The term of highest weight in C^iY) is y\^, ■ ■ ■ y^; that is, 

(1) C,^{Y) = dfcyiS • • • y^ +terms of lower weight, 

where d^ is a constant. 


(a) Ck{Y) is an eigenfunction of the differential operator Ay given by 

( 2 ) 


,2 92 


Z_yi=l tJl Qy'4 ' Ai=l Aj = l, 


J¥=i Vi-Vj 


a 

dyi ■ 


(Hi) As K varies over all partitions of k the zonal polynomial have unit coefficients in the 
expansion of (trY)^; that is, 

(3) {trYf = (yi + • • • + ymf = Ek C^{Y). 


Immediate consequence of Definition [5] is the following important equality 


etr(X) = 


^tr(A:)*^ 

^ k\ 

k=0 


2^2^ k\ 

k=0 K 


(24) 
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Let X,Y e then 


c^{x)Cr{Y)= ^ e;'^c;’^{x,Y) 


(25) 


where 6'^Y = 


C^(Im.) ■ 

For any X,Y G Sm, we have (Gross and Richards, 1987) 


[ C^{XHYH')dH = 

J0(m) 


C^{X)C^{Y) 


I0(m) Ci^{Im) 

For any A G Sm, we have (Gross and Richards, 1987) 


IO<X<Ir, 


det(X)“-i(™+i)a(AX)dX = 

(a+2(m + l))^ 


Let oi,..., Op and bi,... ,bq be complex numbers, such that for 1 < i < p and I < j < q, 
bi > {j — l)/2. Then the hypergeometric function of one matrix argument is defined as 


oFq{(ll , . . . , CLp] bl, . . . , bq] X) — ^ ^ ^ ^ 


k=0 K 


(Q'i)k) • • • ) (®p)k Cf^i^X) 

(^i)k) • • • ) {bq) 


qjK 


k\ 


where denotes the summation over all partition k, k = {ki,, km), /ci > /c 2 > ... > 0, of 
k, and the generalized hypergeometric coefficient (&)«, is given by 

m. . 

(&)«: = n ( ) ’ ib)k = b{b + l)...{b + k-l), (6)o = 1. 


2=1 


The multivariate gamma function which is frequently used alongside is defined as 


F^(a) = / det(X)“-^(”^+^) etr(-X)dX = tt 




2=1 


1 


JJr a--(f-l) , 


where Re(a) > (m — l)/2. 

Multivariate beta function is defined as 

Bm{a, b)= [ det(X)“-i(-+i) _ x)^-5(-+i)dX = 

Jin, r^(a + b) 

where Re(a),Re(6) > {m — l)/2. 


Definition 6 The Laplace transform of the matrix valued function f is given by 


g{Y) = Lf{X) = f etT{-XY)f{X)dX 

Sm 


(26) 


Definition 7 (Press, 1982) A random matrix V G Sm is said to have the non-singular Wishart 
distribution with scale matrix Y G Sm and n degrees of freedom, m < n, if the joint distribution 
of the distinct elements of V is continues with density 

f{V) = cdet(5])-5^det(F)5(”-™-^)etr 

where c~^ = 2~Fm (§)• It is denoted by V ^ WmiY,n). 

Further if we take U = V~^, then U follows the inverted Wishart distribution with scale 
matrix Y and n degrees of freedom denoted by U ^ IWmi'Y,n) with the following density 


--Y-^V 

2 


giU) = cdet(5])5^det([/)-5^-5(™+i)etr 




-1 
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Lemma 15 (Teng et al, 1989) Assume Z is an m x m symmetric matrix, X is an m x m 
complex symmetric matrix with Re(X) G Sm and h is a real function over M"*". Then 



XW)dW 


(a)^r^(a)7fc(a) 

r(am + k) 


det(X)-“C,(ZX-i), 


where Re(a) > (m — l)/2 and 

Ikia) = [ 
Jw.+ 


(27) 


9 Acknowledgements 

We would hereby acknowledge the support of the StatDisT group. We would also want to thank 
Prof. Srivastava (University of Toronto) for his help in Theorem 14. This work is based upon 
research supported by the UP Vice-chancellor’s post-doctoral fellowship programme. 


10 References 

A. Y. Abul-Magd, G. Akemann and P. Vivo, Superstatistical generalizations of WishartLaguerre 
ensembles of random matrices, J. Phys. A: Math. Theo., 42 (2009), 175-207. 

S. Adhikari, Generalized Wishart distribution for probabilistic structural dynamics, Comp. 

Mech., 45 (2010), 495-511. 

G. Aime and R. Soummer, Influence of speckle and Poisson noise on exoplanet detection with a 
coronograph, in EUSIPCO-Of (L. Torres, E. Masgrau, and M. A. Lagunas, eds.), (Vienna, 
Austria), (2004), 509-512. 

T. W. Anderson, An Introduction to Multivariate Statistical Analysis, 3rd Edition, John Wiley, 

New York, (2003). 

T. W. Anderson and K. T. Eang, On the theory of multivariate elliptically contoured distribu¬ 
tions and their applications. Technical Report, No. 54, Department of Statistics, Stanford 
University, Galifornia, (1979). 

M. Arashi, A. K. Md. Ehsanes Saleh, Daya K. Nagar, S. M. M. Tabatabaey and H. Salarzadeh 
Jenatabadi, Bayesian Statistical Inference Eor Laplacian Glass of Matrix Variate Elliptically 
Gontoured Models, Comm. Statist. Theo. Meth., Accepted (2013). 

W. Bryc, Gompound real Wishart and q-Wishart matrices, arXiv.OSOd.^Olfvl [math.PR] (2008). 
Jose A. Dfaz-Garcfa, Ramon Gutierrez-Jaimez, On Wishart distribution: Some extensions, J. 
Lin. Alg. 435 (2011), 1296-1310. 

A. P. Dawid, and S. L. Lauritzen, Hyper-Markov laws in the statistical analysis of decomposable 
graphical models. Ann. Statist. 21 (1993), 1272-317. 

J. Pan, Non-central Gochran’s theorem for elliptically contoured distributions, Acta Math. 

Sinica, 3(2) (1986), 185-198. 

G. Erahm, Ceneralized Elliptical Distributions: Theory and Applications, PhD Dissertation, 
Universitat zu Koln, (2004). 

Prancisco J. Garo-Lopera, Graciela Gonzalez-Parias, N. Balakrishnan, On generalized Wishart 
distributions - I: Likelihood ratio test for homogeneity of covariance matrices, Sankhya A, 
DOI 10.1007/sl3171-013-0047-7, (2014). 

E. D. Eeigelson and G. J. Babu, Modern Statistical Methods for Astronomy With R Applications, 
Gambridge, UK, (2012). 

A. K. Gupta and D. K. Nagar, Matrix variate distributions Ghapman &: Hall/GRG, (2000). 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, 7th Ed., Academic 
Press, USA, (2007). 

K. I. Gross and D. ST. P. Richards, Special functions of matrix argument I: Algebraic induction 

zonal polynomials and hyper geometric functions, Trans. Amer. Math. Soc., 301 (1987), 
475501. 


21 



Anis Iranmanesh, M. Arashi, D. K. Nagar and S. M. M. Tabatabaey, On Inverted Matrix Variate 
Gamma Distribution, Comm. Statist. Theo. Meth., 42(1) (2013), 28-41. 

G. Letac and H. Massam, The normal quasi-Wishart distribution. In: Viana MAG, Richards 

DStP (eds) Algebraic Methods in Statistics and Probability. AMS Contemp Math 287 
(2001), 231239. 

E. Lukacs, and R. G. Laha, Applications of Characteristic Functions, Griffins Statistical Mono¬ 
graphs & Courses, 14. Charles Griffin &: Co., Ltd., London, (1964). 

Rob J. Muirhead, Aspect of Multivariate Statistical Theory, 2nd Ed., John Wiley, New York, 
(2005). 

S. Munilla and R. J. C. Cantet, Bayesian conjugate analysis using a generalized inverted Wishart 
distribution accounts for differential uncertainty among the genetic parameters an appli¬ 
cation to the maternal animal model, J. Anim. Breed. Genet., 129 (2012), 173187. 

S. J. Press, Applied Multivariate Analysis, Holt, Rinehart &: Winston, New York, (1982). 

A. Roverato, Hyper-Inverse Wishart Distribution for Non- decomposable Graphs and its Appli¬ 

cation to Bayesian Inference for Gaussian Graphical Models, Scandinavian J. Statist., 29 
(2002), 391-411. 

M. S. Srivastava and C. G. Khatri, An Introduction to Multivariate Analysis, North-Holland, 
Amsterdam, (1979). 

B. C. Sutradhar, and M. M. Ali, A generalization of the Wishart distribution for the elliptical 

model and its moments for the multivariate t model, J. Mult. Anal, 29 (1989), 155-162. 

C. Teng, H. Fang, W. Deng, The generalized noncentral Wishart distribution, L. Math. Res. 

Exp., 9(4) (1989), 479-488. 

J. Y. Tourneret, A. Ferrari, G. Letac, The noncentral Wishart distribution; Properties and 
application to speckle imaging, in: Proc. IEEE Workshop on Stat. Signal Proc., (2005), 
18561860 

H. Wang and M. West, Bayesian analysis of matrix normal graphical models, Biometrika, 96(4) 

(2009), 821-834. 

C. S. Withers and S. Nadarajah, log det A = tr log A, Int. J. Math. Education Sci. Tech., 
41(8) (2010), 1121-1124 

C. S. Wong, T. Wang, Laplace-Wishart distributions and Cochran theorems, Sankhya 57 (1995), 
342-359. 


22 



